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ABSTRACT 

Global stability robustness with respect to unmodeled dynamics, arbitrary 
bounded internal noise, as well as external disturbance is shown to exist for 
a class of discrete-time adaptive control systems when the regressor vectors 
of these systems are persistently exciting. Although fast adaptation is 
definitely undesirable, so far as attaining the greatest amount of global 
stability robustness is concerned, slow adaptation is shown to be not 
necessarily beneficial. The entire analysis in this paper holds for systems 
with slowly varying return difference matrices; the plants in these systems 
need not be slowly varying. 
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1 . INTRODUCTION 


Persistency of excitation (PE) has been claimed to be essential for guar- 
anteeing some robustness in adaptive control (AC) systems ever since it was 
proven to ensure robust exponentially fast identification [1,2]. But, to 
date, no formal proof exists that PE guarantees global stability robustness 
for discrete time systems in general. In [3] it is shown that if the refer- 
ence input is PE and the signal-to-disturbance ratio is large, then model 
reference adaptive control (MRAC) system signals will be globally bounded. In 
[4] continuous time model reference adaptive control systems are shown to be 
globally stable with unmodeled dynamics if the reference model is strictly 
positive real and the regressor vector is persistently exciting. Attempting 
to extend the same to discrete time MRAC systems failed due to the fact that 
discrete time MRAC systems do not exactly parallel continuous time MRAC 
systems [5, Sec. 5.3]. As a result, only local stability is proved and slow 
adaptation has to be assumed. No direct and explicit proof has yet been shown 
that persistency of excitation guarantees some global stability robustness 
against nonzero unmodeled dynamics and bounded internal as well as external 
(with respect to the plant) noises for an entire class of discrete time AC 
systems. 

Perhaps fueled by the above bold claim together with a lack of any formal 
proof that PE can guarantee robustness in general, many researchers resorted 
to showing examples of problems of adaptive control systems and developing ad 
hoc remedies for these problems [6,7,8]. It is well known in (nonadaptive) 
robust controls research that a controller design based on the nominal plant 
properties has only limited robustness against unmodeled dynamics in the plant 
[9]. Showing robustness problems of AC systems demonstrates that the limit 
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extends to AC systems as well, and one should not expect any AC system to 
perform satisfactorily under all unpredictable situations. The remedies in 
[6,7,8] have value in solving the problems specifically exhibited provided 
their required prior informations (e.g., solution starts within certain 
bounded set [6], etc.) are verified. However, if an AC system does possess 
some inherent global stability robustness, due to say PE, whether further 
incorporating these ad hoc remedies would increase or decrease such basic 
robustness is a question that remains to be answered. 

The contribution of this paper is in proving directly and explicitly that 
when the regressor vector of a discrete time AC system of a particular class 
(maybe with a fast time-varying plant, as long as the system return-difference 
matrix [9] is uniformly bounded and sufficiently slowly varying) is PE, the 
system global stability is inherently robust against nonzero unraodeled 
dynamics expressed in terras of stable factor perturbations [9] as well as any 
bounded noise internal or external to the plant. It is also shown that if the 
adaptation gain approaches zero (slow adaptation) the margin of such robust- 
ness approaches a constant. However, this does not imply that slow adaptation 
necessarily improves global stability robustness. This is justified by 
analyzing a special case when the regressor vector is a scalar. The conclu- 
sion obtained based on this analysis is that, while increasing the adaptation 
rate indefinitely reduces the robustness for sure, the maximum robustness of 
the AC system may be achieved at a finite adaptation rate. 

Thus, if the popularity of using the averaging technique in AC system 
stability analysis and the abundance of results obtained based on slow adapta- 
tion, e.g., [5,10], in any way imply a strictly non-negative relationship 
between slow adaptation and global stability robustness, this result shows the 


contrary. 



This paper is organized as follows. In the following section, some nota- 
tions and useful facts on robust control and linear time-varying systems are 
first listed for ease of referencing later. In Section 3, the class of adapt- 
ive control systems under consideration is developed to the extent relevant to 
the derivation of the results of this paper. A general error model is then 
developed for this class of AC systems in Section A, followed by the inherent 
global stability robustness proofs assuming persistency of excitation of the 
regressor vector. In Section 6, the effects of the adaptation gain on the 
inherent robustness are analyzed. Finally the results are summarized and the 
paper is concluded with a discussion. 


2. PRELIMINARIES 

Let z denote the transfer function of the unit-delay operator. Let 
R[z] denote the ring of polynomials in z with real coefficients and R(z) 
denote the field of rational functions assoicated with R[z]. Let Z e R(z) 
denote the set of causal stable rational functions of z, i.e., functions with 
poles strictly outside the closed unit disk. Let M(Z) denote the set of 
matrices whose elements belong to Z. 

The plants considered in this paper are lumped linear causal discrete- 
time MIMO systems whose transfer matrices (or sequences of transfer matrices 
if systems are time-varying) belong to M(R(z)), where M denotes matrices. 

Definition 2.1: Let U(Z) denote the set of unimodular matrices in 


M(Z). It is defined by 
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U(Z) = {F e M(Z) : F 1 e M(Z)} . 


Fact 2.2 : M(R[z] ) C M(Z) . 

Fact 2»3 : Let E, F e M(R[z]). If E, F are left or right coprime 

in M(R[ z ] ) , then they are also left or right coprime in M(Z) respectively. 

Let a plant P e M(R(z)) and a controller C e M(R(z)) be connected 
as shown in Figure 1 . 



e 
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Figure 1. Standard feedback control system 

Let (N,D) e M(Z), ND * = P, be a right coprime (in M(Z)) factorization 
(r.c.f.) of P and (X,Y) e M(Z), Y -1 X = C, be a left coprime (in M(Z)) 
factorization (l.c.f.) of C [9]. 

Fact 2.4 : [9] C stabilizes P in the bounded-input bounded-output 

(BIBO) sense if and only if 


XN + YD e U(Z). 
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XN + YD is referred to as the return-difference matrix of the system in 
Figure 1. BIBO stability is often referred to as internal stability also. 

Fact 2.5: The poles of the system in Figure 1 are the zeros of XN + YD. 


Let E e M(Z) be expressed as the infinite sum 


E(z) = l E z 1 . 
i=0 


Definition 2.6: Define the norm II • II _ : M(Z) + R by 

3 


i eii = y he. n. 
a j=o J 10> 


where H • H ^ denotes the induced infinity norm [11]. 


Fact 2.7 : Let E e M(Z). Let y^fE} denote the i^-gain of E 

[21, Section 3.2]. Then, for the systems under consideration, 


yjE} - »En a . 


Fact 2.8: [12] Let U e U(Z) and V e M(Z). If II V - Ull „ < 1/llU 1 | - 

" — 1 3 3 

then V e U(Z). 


Fact 2.9: [25,19] Let E, F, G e M(R[z]). If E, F are right coprime 


in M(R[z]), then there exists X, Ye M(R[z]) 


such that 
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XE + YF = G. 


3. THE ADAPTIVE SYSTEM 

The class of AC systems considered in this paper consists of MIMO AC 
systems with direct signal-feedback connections (as opposed to only parameter- 
feedback via the parameter adjustment process) from the outputs through the 
controllers to the inputs of the plants. The primary objective of adaptation 
for these systems is to satisfy asymptotically the necessary and sufficient 
condition for internal stability stated in Fact 2.4. An example of an AC 
system with only parameter-feedback (no explicit signal-feedback) appears in 
[13]. A general block diagram of this class of AC systems is shown in terms 
of an r.c.f. (N,D) of the strictly causal plant P and an l.c.f. (X,Y) of the 
controller C in Figure 2, where d denotes disturbances 


d n 



Figure 2. General AC system 

external to the plant and n noises internal to the plant. This class of 
systems includes those studied in [14,15,16,17,12]. A general AC methodology 






t 
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for achieving the asymptotic internal stability objective can be developed 
roughly as follows. Complete details can be found in [14] for a system in the 
less general observer-controller structure and in the unit-advance operator 
representation, and in [16] for a system in exactly the same structure as that 
in Figure 2 and in the unit-delay operator (z) representation. 

Ignoring n, d, the identifier, the time index k and the estimation 
symbol consider the feedback system in Figure 2. Suppose P e M(R(z)) 

is not precisely known and the objective is to find a C e M(R(z)) such that 
the feedback system is BIBO stable with poles at the zeroes of 

U x = diagfu^ e U(Z), e R[z] ¥ i. 

Assume N,D e M(R[z]) and are coprime in M(R[z]) from now on. By Facts 
2.2, 2.3, 2.4, and 2.5, if one can find a pair (X,Y) e M(R[z]) such that 

XN + YD = U , (1) 

then implementing C by Y will achieve the above stability objective. 

To find (X,Y) adaptively through identification, an error function e 
needs to be defined. The convergence of this function to zero must result in 
a solution pair for (1) to be identified. To find such e, let x := D ^u so 
that 

y=ND^u=Nx, u = Dx . 


Multiplying (1) by x results in 



Since x cannot be measured directly, it must be estimated also. To do this, 


note that since N and D are right coprime in M(R[z]), there exist A^ , 

Bj e M (R[ z ] ) (to be called auxilliary-controller) such that 

(3) 

(4) 

Now, multiply (4) by x to get 


and, if (1) holds 


AjN 


+ BjD = I 


XN + YD 


U 1 (Aj^N + B l D) . 


Xy + Yu = Uj(A^y + B^u). 


(5) 


Finally, define e by 


e = UjCAjy + Bju) - (Xy + Yu). (6) 

But, this is not the end to completely defining e. This is because if we 
apply a linear estimation algorithm to make e(k) converge to zero now, the 
solution (A^ , Bj) can take a form not necessarily satisfying (3). 
Fortunately, it can be shown [16] that by imposing degree constraints on X, 
Y, Aj, Bj, Uj» and restricting the constant term of b 1 to be lower- 
triangular with l's in the diagonal, the (A^, Bj) that solves (6) can be 
forced to take values such that 
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|A l N + B j^D | = 1 

when e (k) converges to zero. Therefore, although the to-be-identified 
solution pair (A^ , Bj ) may still not satisfy (3), it will not change the final 
feedback system poles from the desired values which are the zeros of Up 

Note that if P were stable to begin with [12] , D would be unimodular 
in U(Z), and (1) can be equated to D instead of Up i.e., 

XN + YD = D. (7) 

In this case, multiplying (7) by x yields 

Xy + Yu = u. 


Defining e by 

e = u - (Xy + Yu) (8) 

a much simpler error function for identifying a controller for the system is 
obtained. The final poles of such a feedback system would be the same as the 
poles of the plant. 

If instead P were stably invertible with known right interactor 
matrix J = z -0 H, H a causal polynomial matrix [16], then JN e U(Z) and 
(1) can be equated to JN instead of Up i.e., 

XN + YD = JN = z~ 5 HN. 


(9) 
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Multiplying (9) by z^x = z^D yields 

z^ Xy + z^Yu = Hy. 

Defining e by 

e = Hy - (Xz 6 y + Yz 6 u) (10) 

again results in a linear function in terms of the unknowns (X, Y) which can 
be identified by applying any linear estimation algorithm. If the identi- 
fied X and Y actually satisfy (9) the feedback system will have all poles 
at z = «. If the system is further implemented as shown in Figure 3 instead 
of 
v 

y 


Figure 3. Delayed-tracking system 

as in Figure 2, the delayed-tracking scheme of [17] results. 

For the rest of this paper, we will consider only the general case when 
the plant is not necessarily stable or stably invertible and the error func- 
tion is defined by (6). 

To identify the unknowns Aj , Bj, X, and Y using (6), note that since 
is diagonal the unknowns can be identified row by row independently. 
Hence as far as identification is concerned, the MIMO system being considered 
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ls no different from a SISO system. For integrity and generality, however, we 
shall keep the MIMO nature of the system. Let a^, bp xp yp 
i = 1, •••, n denote the rows of Ap Bp X, Y respectively. Let 



U 


U l X 

V 


u I 
n 


Then (6) can be rewritten as 


< 

r-H 

II 

U) 

U 1 B 1 

X Y] [y' u' 

-y' -u']' 

= [AU 

BU 

X Y] [y' u' 

~r -u']' 

= [A B 

X Y] 

[(U y y (U u)' 

1 

\ 

1 

c 

\ 


Let 0^ denote the ordered vector of parameters of [a^ b^ y^] to be 

estimated. (Note that usually some of the parameters of [a^ b^ x^ y^] 
are assumed known and need not be estimated.) Let 0 ^(k) denote the k-th 

estimate of 0^ and <t>^(k) denote the k-th ordered regressor vector 

A 

corresponding to 0^(k) and made up by the measurements of ~[u^y' u^u' 
-y' -u ' ] ' . Furthermore, let 
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1 

<D 


QjOO 

0 2 

a 

0 9 (k) 

• 

, 9(k) : = 

L. 

m 

• 

• 


• 

1 

c 

<D 

J 


6 (k) 

L n J 


r (k) 


r 1 (k) 

r 2 (k) 


o 


4> (k) 


I j ( k ) 

<J > 2 (k) 


A (k) 


X (k) 

X 2 (k) 


k j 


0 




(J) 3 (k) J 

0 1 


( k ) 

n 


where T^(k)'s are symmetric positive definite gain matrices and X i (k)^s 
are positive scalars. The equation error identification algorithm to be used 
is assumed to be in the form 


e(k+l) = e(k) + r(k) 4 >(k) (A(k) + <T(k) r(k) 4. (k) ) X e(k), 


e(k) = i|i(k) - 4'(k) 9 (k) . 


( 12 ) 


Here iji(k) denotes the vector of k-th values of signals generated by the 
not-estimated parameters of [A B X Y] operating on their corresponding signals 
according to (11). Identification algorithms having precisely this form 
include at least the recursive least squares and the normalized gradient algo- 
rithms [5], Other algorithms having different forms may also be used. 
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In case the plant were stably invertible and the simpler error function 
(10) were used for identifying the controller, the identification algorithm 
would be assumed to be in the 6-interlaced form 


d (k+1 ' 


0 (k-<5 ) 


r(k-6) 4>(k-5)(A(k-6) + <J>'(k-6) r (k -5 )$ (k -<5 )) 

e(k-d) = f (k— 6 ) - <(.'(k-6) 0(k-5). 


L e ( k— 6 ) 

(13) 


Note that this is different from the algorithm used in [5, Section 5.3.2] in 
that e (k-6 ) is a function of 4>" (k— 6 ) 0(k-6) instead of <|>'(k-6) 0(k). 
This is important because the signal values in <f> (k— 6 ) are exactly what 

the k-6th estimate (x , y ) of the controller operates on when it 
is implemented as the controller for one period of time. This is also the key 
factor in eliminating the s low-adaptation requirement in proving global 
stability. 

(13) is equivalent to 


0(k+c) = 0 (k) + r(kH(k)(A(k) + <{>'(k)r(k)(J>(k)) -1 e(k) , (14) 

e(k) = *|>(k) _ <(»'(k)0(k), ? 1 . 

Since this is a more general form of identification algorithms than (12), we 

shall assume that (14) will be used rather than (12) from now on. 

The controller C will be assumed to be updated each time (X,Y) is esti- 

— 1 A k 

mated and implemented as (Y ) X on line. Hence the implemented AC 
system will have exactly the structure shown in Figure 2. We are now ready to 
derive a complete error model for the class of AC systems under consideration. 
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4. THE COMPLETE ERROR MODEL 

Suppose P has been undermodeled so that the pair (N,D) is only a 
nominal (or tuned [4]) right coprime (in M(R[z])) factorization of P. 
Correspondingly, (X,Y) and (A,B) are only left coprime (in M(R[z])) factori- 
zations of the nominal controller and nominal auxiliary controller of P 
respectively. Let (N^, D*) e M(R[z]) be a true r.c.f. (coprime in 
M(R[ z ] ) ) of P. That is, N *D* * = p. Then the stable factor perturba- 
tions (N^ , D^) , 

N* - N - N*, D* = D - D*, (15) 

are obviously due to the unmodeled dynamics in P. The purpose of this paper, 
as mentioned in the introduction, is to prove that if the adaptive methodology 
developed in the preceding section for the nominal plant ND * is applied to 
the true plant as shown in Figure 4, the resulting AC system will 
stay globally BIBO stable if II [N', is sufficiently small, r, d, n 
are uniformly bounded, and the regressor matrix 4> (k.) is PE. This will be 
done through error model analysis. 


d n 



Figure 4. True AC system 
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To develop an error model for the complete undermodeled AC system, first 
subtract 9 from (14) and let 0(k) := 0(k) - 0 to yield 


0(k+U =0(k) + r(k)<|>(k) (A(k) + <J>'(k)r(k)<(>(k)) ^(k). 


e (k) = i(;(k) - <{> ^ (k) 0 (k) • 


(16) 


Clearly, 0 (k) must be the ordered vector of coefficients of 


* k * It A k ~k ~k ~k ~k 

[A -A B -B X -X Y -Y] = : [k B* X* Y K ] . 


A k 

corresponding to those of [A 

pointing out that coefficients of 

* k A k 

estimated parameters of [A B 
tion. Now note from (14) that 


B k X k Y k ] estimated • Here it is worth 

[A k B k X k Y k ] corresponding to the not- 
A k A k 

X Y ] are all zero due to the subtrac- 


Mk) “ $"(k) 0(k+£) * i|;(k) - (j>"(k)0(k) 

- <r(k)r(k)4>(k)(A(k) + <|)'(k)r(k)<{.(k)) _1 (iKk) -♦'(k)e(k)) 
= A(k)(A(k) + <t>'(k)r(k)<j>(k)) -1 (iKk) - $'(k)0(k)) ; 


therefore 


0(k+c) = 0(k) - r(k)<Kk)A 1 (k)(*'(k)e(k+c) - *(k)) 


(17) 
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e (k+c ) — -r(kH(k)A l (k)(<rOOe(k+?) -*00). (18) 

1 - z Q 

Let (A^, B*, , B^, X A , Y fc ) e M(R[z]), A*, B^ in exactly the same 

block forms as A, B respectively, be a solution of 

X * N * + Y * D * " A * UN * + B * UD * 

= Vl* N * + U^D* (19) 

=: V* e U(Z) 

closest to (A, B, Ai , Bi , X, Y) in the H • II norm. Note that a solution 
of (19) always exists because N*, D* are right coprime in M(R[z]) (Fact 

2.9). By Fact 2.4, this implies that (X^, Y^) and (A B^) are a true 

controller and a true auxiliary controller of P respectively. Redefine x 
by D A *u then (referring to Figure 4) 

y = N^(x + n) , u = D^x. 

Multiplying (19) by x yields 

X*(y - N*n) + Y^u = A*U(y - N*n) + B*Uu = V*x. (20) 


Let 



then by the various definitions 
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t|>(k) - 4> ' ( k.)0 (k+£ ) = A k+C Uy(k) + B k+? Uu(k) - X k+C y(k) - Y k+C u(k). 

= (A* + A k4? ) Uy(k) + (B*+ B k+C )Uu(k) 

- (x* + x k4? )y(k) - (Y* + Y k+? )u(k) (21) 

- (X* - A*U)N*n(k). 


Define V and by 

V- XN* + YD* e M(Z), (22) 

V a = AUN* + BUD* e M(Z). (23) 

Remark 1 : Note that if the plant is not undermodeled X, Y, A, and B 

exist such that V = V g and the global stability issue can be much more 
easily addressed. With undermodeling, V and can only be close to each 

other (in terras of small II V - V II ) in general if the degree of under- 
modeling is slight. 

Multiplying (22) and (23) by x, one gets 

Vx = X(y - N*n) + Yu, (24) 


V a x - AU(y - N*n) + BUu 


(25) 
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Subtracting (20) from (24) and (25) 

X*(y - N*n) + Y*u = (V - V*)x =: Vx, 

A*U(y - N*n) + \uu = (V - V*)x =: V x. 

3 . 3 

Substituting into (21) 

t|> (k) - 4> " ( 1c )9 (k +5 ) = A k+? Uy(k) + B k+? Uu(k) + AUN.n(k) + V x(k) 

* 3 

-X k+C y(k) - Y k+? u(k) - XN*n(k) - ^x(k) (26) 

= -<{.'(k) 0 (k+c) - (V - V )x(k) - (X - AU)N.n(k). 

3 * 

Now, from Figure 4 it is obvious th3t 

[X k Y k ] [y' u']' = [X k I] [r' d'V . 


Therefore , 

Vx = [X Y][y' n'V - XN*n 

= -[X k Y k ][y' n'V + X k r + d - XN*n. 

Let £(k) denote a matrix entirely similar to <|>(k) except made up of 

the measurements of r and d instead of those of y and u respec- 


tively. Let also 
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r ■ 

O 


' 0 

II 

sT 

0 

I 

T • = 

’ L l • 

0 

I 


I- 


0. 


Then 


= V 1 ([A k B k X k Y k ] (T L 


r _ 


r-u . yl 

~r 



-d 

- M, 

-u .u 


i 

r 

1 


. d. 


1 


+ X(r - N.n) + d), 


x(k) = V ^'(lOT - <T(k)M)0(k) + X(r(k) - N*n(k)) + d(k)) 


(27) 


where T,M are diagonal matrices corresponding to Tj, Mj and mask out the 
elements of 0 (k) corresponding to the parameters of A k B k , Y k 
appropriately. Thus, a complete error model is finally developed. It is 
described by (18), (26), (27) and is shown in Figure 5. 




















- 21 - 


5. ROBUST GLOBAL STABILITY 

The global stability proof of the error system in Figure 5 will be 
analyzed by simply applying the well-known small gain theorem and assuming 
persistently exciting <J> ; that is, there exists a, 3, N, N J> 1, 

0 < a < 6 < 09 such that for all k > (N-l)s 

N-l 

01 > l <t>(k-ic)<T(k+ic) > a I. (28) 

i=0 

The error system will be considered as consisting of three subsystems 

h I* 0(k-G), 

S 2 : 9(k-c) |+ (a,b) , 

S 3 : c |+ g. 

We'll begin with the well-posedness issue. Referring to Figure 5, Sj 
relates h to 6 (k-£ ) as follows: 

0 (k+o = — - r(kH(kH _1 (k)(h(k) + ♦'(k)e(k+c», 

l-z’ 

e(k+o - - (i + _i_ r(kH(k)A _1 (k)<|.'(k))" 1 r(k)<|>(k)A" l (k)h(k) 

1-z 5 1-Z g 

= - (I - z C I + r(k)<|>(k)A~ 1 (k)()>'(k)) _ 1 r(k)<t»(k)A~ 1 (k)h(k). 

0(k+O = (l+r(k)<t>(k)A _1 (k)(j)'(k)) -1 (9(k) - r(k)<j>(k)A -1 (k)h(k)) 
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I 


= r^k)(i + r 1 ^Ck)4> (k)A -1 (k)<t> (k)r^k)) 1 r" L/ ^k)(e(k) - r(k)<Kk)A -1 (k)h(k)) . 
S2 and S3 can be described by the relations 

(b - a) (k) = (£'(k)T - <T(k)M)z ? e(k+5)> ( 29 ) 

g(k) = (V - V fl ) V _1 c(k). 

Signals h and c are related to d, f, and p (which are uniformly bounded 
if r, n, and d are uniformly bounded) as follows: 

c = d + f + S^S^(p + S^c) 

c = (I - s 2 s 1 s 3 )" l (d + f + s 2 s lP ) 

h = (p + S 3 c) = p + S 3 (I - S 2 S 1 S 3 )~ 1 (d + f + S 2 S lP ). 

Since A(k) > 0 , I + (k)(J> (k)A 1 (k)<t» " (k)r (k) >_ I, therefore 

(i + I^ 2 (k)<f>(k)A 1 (k)<|> '(k)r ^ 2 (k)) 1 is well-defined and is well-posed 

as long as 9 ( 0 ) is bounded. S2 is well-posed as long as <(> (k ) and 

5(k) (hence r and d) are uniformly bounded. is well-posed because 

V, V e M(Z) and the assumption that P is strictly causal (hence N* is 
a 

strictly causal) imply that the constant term of V = XN* + YD* is non- 
singular and so V * e M(R(z)). Therefore the individual subsystems are 

well-posed. For the entire system, since S2 consists of at least one unit 

5 > 1, (I - S 2 SjS 3 ) * is causal and well-defined. 


of delay due to 
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Therefore, if d, f, p e l then k, c e i and the small gain theorem 

00 0 CO 0 

can be applied to the error system. 

The homogeneous part of is made up by the relation 

8(k+c) = (I + r(kH(k)A -1 (k)«|.'(k)) _1 0(k). (31) 

This relation has been shown to be exponentially stable for many linear esti- 
mation algorithms [18,19] in which T(k) is updated in certain ways and the 
elements of A(k) are constant and initialized between (0,1] provided that 
the PE condition (28) is satisfied. Assuming one of these algorithms is used 
and again (28) holds, then the exponential stability of (31) implies that 
is BIBO stable [20, Lem. 2.2] and < ». 

With Sj being BIBO stable the following is true. 

Lemma 2: Provided (28) holds and r, d are both uniformly bounded by 

say, §2 is BIBO stable. 

Proof : (28) implies that 

"<J>(k)I ito < 8 1/2 < - V k; 

that is, <|>"(k)M is BIBO stable. Since £(k) is made up by the measure- 
ments of the uniformly bounded r and d, £'(k)T is also BIBO stable. 

Therefore, S 2 is BIBO stable. 
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Hence, S2SJ is BIBO stable with 

Y »{ S 2 S lJ - ( max {^ 2 > =: 0 4 6 i < 08 * 

By the small gain theorem [11, Theorem 3.2.1], the entire error system 
will be BIBO stable if Y 00 {S 3 } can be made smaller than l/y^fs^}. 
This is always possible due to the following theorem. 

Theorem 3: yj S3} = II (V - V a )V -1 H g is smaller than l/yjs^} if 

ll[N* D*]'ll a is sufficiently small. 

Proof : By definition, 

V = XN* + YD* - XN + YD - (XN* + YD*) e M(Z) 
where XN + YD e U(Z). Therefore, 

XN* + YD* = XN + YD - V. 

By Fact 2.8, if II [X Y] [n; 'fi']*’! < l/« (XN + YD) -1 n & then V e U(Z), 

V * e M(Z) and 

n v -1 !! =: 0. < ~. 

a j 


By definition again, 


V 



V - V 

a 
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= XN* + YD* - (AUN* + BUD*) 

= [x-au y-bu] [n; D']'. 

Applying (4) and noting that AU = U^A^, BU = U^B^ yields 

v - v a = -[x-au y-bu][n; d;]', 

II V - V II < II [X-AU Y-BU ] II II [NT Dll'll • 
a a — a * * a 

Since X, A, U, Y, B e M(Z), II [X-AU Y-BU ] II =: 3 , < ». Therefore, 

a o 

b (v - v )v _1 ii < ii v - v n h v _ 1 n 

a a — a a a 

< b 5 b 6 iin; D;n a . 

Hence if II [ N' 1 ^ n a is smaller than 

rain{ 1/ (II (XN + YD) _1 II a ll[X Y] n a ) , l/CU^yJS^} )} 

then II ( V -V )V _1 D < 1/y (S„S,). 

a a 00 z I 


Thus the class of AC systems under consideration is proved to be robustly 
globally stable against unmodeled dynamics and bounded n and d using the 
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most conservative small gain theorem and assuming persistency of excitation of 
the regressor vector. It is expected that the margin of such robustness will 
be less-conservative in practice than that the proof of Theorem 3 has perhaps 
impressed upon the reader. 

Remark 4 : Note that the key to the existence of a nonzero global 
stability robustness margin is the finiteness of y^CS^S^). The PE condi- 
tion only ensures this; it is not a necessary condition for global stability 
robustness. 

The comment preceding Theorem 3 implies that the smaller Y oo{ S 2 S l^ is 
the more robust the global stability of the AC systems under consideration 
will be. In the next section, we offer a qualitative analysis of how the 
degree of persistency of excitation, i.e., the values of a, 6 (28), and the 

gains T(k), A(k) affect the smallness of y { S t } . 

°° /. I 


6. ROBUSTNESS VS. DEGREE OF PE AND ADAPTATION GAINS 

Due to the difficulties of analyzing the effects of time-varying matrix 
gains on the robustness in general, we will consider only SISO AC systems and 
assume T(k), A(k) to be equal to the scalar constants y, X respec- 


tively in this section. 
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Recall from Section 5 that Sj relates h to 0(k-t;) by (28) or 


9(k+c) = (i + y 4,(k)r(k)) _1 (0(k) - I <j>(k)h(k)) 


- (I + y 4.(kH'(k)) -1 0(k) - (i + y ♦(k)*'(k))" 1 2- ♦(k)h(k) 


(33) 


m 


( n (I + 1 *(k-icH'(k-ic)) ^©(k-mc) 
i=0 A 


m J 


I ( n (I + f (Kk-ic)^'(k-ic)) 1 ) f ♦(k-jc)h(k-jc) 
j=0 i=0 A 


where m := fk/c"). Let t := [j/Nj, since (i + <j>(k-i£)<f>'(k-ic))~* < I 

Vk > i c 

I0(k+c)ll < Il0(k-mc)il + 


(34) 


m t AN-1 v 1/ 

+ l ( n II n (I + I ^(k-ic)^'(k-ic)) »,J f B /2 |h(k-jc)| , 

j=0 A=1 i=(A~l)N 


where if t < 1 the term within the parenthesis is defined to be equal to 
1. Now, 


v -1 U ' 1)N 

n (i + f <Kk-is)<r(k-is)) - ( n (i + y ♦(k-icH'(k-ic))) 

A J A IT 1 A 


AN-1 

n 

i=(A-l)N 


-1 


i-AN-1 


AN-1 

l 

i=(A-l)N 


- (I + I ♦(k-ic)*'(k-ic) + ••• + (y) N n ♦(k-ic)*'(k-ic))“ 1 . 


AN-1 

n 

i=(A-l)N 


(35) 
-1 


If y is small and approaches 0, then 


AN-1 AN-1 

n (I + I <Kk-ic)<j.'(k-iO ) -1 + (I + 1 l <Kk-ic)<r(k-is)r l 


i«(A-l)N 


i=(A-l)N 
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< (I + f a I) 1 = r-4 I, 

— X X + ya 


Il0(k+?)ll oo < 110 (k-m 5 >11^ + ( l ( — — ) t f8 1/2 ) max |h(k-j?) 


J-0 


j=0,m 


< 110 (k-mi; >11^ + (N B^ 2 l ( - * — ) max |h(k-jc) 

j=0 T j=0 ,m 


10 (k-m£ )ll m + — — — ^ max |h(k-jc)| 

j=0,m 


* lie (k-m? )H w + ^— max |h(k-jc)|. 

j=0,m 


NR i •• 

This implies that if y -*■ 0. From Figure 5 it is obvious 

'=0 < h 


-V 2 

that if llr(k)ll < 0 7 V k then 


Yjs 2 } = Yj(r(k)T - <(.'(k)M)z ? } < k L 8 7 + k 2 e 1/2 


where denote the ranks of T, M respectively. This leads to 

Yj S 2 S il < ^ 0^7 + k 2 ^ 2) - (36) 

Therefore, if y^fs^} < a/(NB^ 2 (k^B^ + k^^ 2 )) =: p then the error system, 
hence the AC system, will be globally BIBO stable for y + 0. In this sense 
the value p may be termed the inherent (due to PE) robustness margin of 
the AC system when y 0. 

The above result as well as the large number of AC system robustness 
results obtained based on slow adaptation (or the averaging analysis) [5] does 
not imply that slow adaptation necessarily increases or decreases the inherent 
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robustness of an AC system. This can be enlightened by analyzing the robust- 
ness margin of a S1S0 system with only one parameter to be estimated. 

Consider (35); if <(> (k.) is a scalar, then all terms under the 
inversion sign are aonnegative. Therefore, if (28) is satisfied. 


JtN-1 2 

n (1 + |-4>(k-ie)<r(k-iO) < (1 +P- + I T A(k)) 

i»U-l)N A A X Z 

2 

y 

where A(k) denotes the nonnegative contribution by the terms in (35) of 

X 2 

orders or higher. Assuming A(k) > n > 0 for all but finite number 

X 

of k, then it is reasonable that a A(n) > 0 exists such that on the 
average 


JIN-1 2 

n (l + \ <Kk-iO<r(k-iO) < (1 + P- + \ A(n))" 1 . (37) 

i=U-l)N A K 


Substituting (37) into (34) and go through similar manipulations as for the 
Y 0 case results in 


119 (k+C )n oo < H8(k-m?)ll oo 




(y + 


aX + YA(n) 


) max 
j=0,m 


|h(k-jc)| 


and therefore 


y { s 1 = W-Tfv + * •) = m % f X. 

]7 X ^ aX +YA(n)^ ^X 


a + A(n) 


(38) 


From (38) it is obvious that 



as y 0, 
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yJh) 


as y 


l/ 2 

But — — is not necessarily the smallest y { S , } . This can be verified 

a “ 1 1 J 

by differentiating (38) with respect to y and equating to zero to yield 


2 a A (n ) + (^-) 2 A 2 (n) + (a 2 - A(n)) - 0. 


(39) 


This equation has only one valid solution at 

L = + ^A(n) 

x A(n) "• Y * 

2 

if A(n) :> a . At this solution point 

N8 ^ 2 

which is smaller than except when A(q) = a . This implies that 

2 

provided A(n) > ot , reaches its lowest value (hence the AC system 

is the most robust against unmodeled dynamics and disturbances) when the true 

v 

rate of adaptation ^ is not the slowest possible. 


7. CASE OF TIME-VARYING PLANT 

Suppose the plant P had been not only undermodeled but also time- 
varying. The development in Section 4 would go through entirely while the 
results of Sections 5 and 6 would still hold true almost completely if one 
simply adds a superscript k to the symbols 
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N*, D*, N, D, X*, Y*, X, Y, A*, B a , A, B, A^, B^, 

Aj, N*, D*, X*. B*, 9, V*, V, V, V a , and V & . 

The results of Sections 5 and 6 will become completely true if one further 

k -1 

assumes that the now time-varying subsystem (V ) (see Figure 5) is 

k -1 

sufficiently slowly varying while satisfying (V ) e M(Z) V k and 
IIN*II , IID k ll , BX k H , IIA k || all uniformly bounded over k. 

3 3 3 . 3 

These are required so that (V k )~^ is BIBO stable [28,29] with finite 
Y w {(V k ) and f, p, <|> are uniformly bounded as assumed. Then 

{v k - V k } can be greater than zero and still satisfy the unity-loop-gain 
condition of the small gain theorem. Thus the global BIBO stability of the AC 
system with a time-varying plant is also proved to be robust. Note that 
although assuming the plant varies slowly would make the assumption that 
(V k ) _1 i S sufficiently slowly varying more easily satisfied, there is no 
apparent need to do so, not even for the sake of achieving the desired per- 
formance as indicated by our preliminary analysis (not included). 


8. RESULTS AND DISCUSSION 

It has been shown that, provided the regressor vector is persistently 
exciting, the class of adaptive control systems with possibly time varying 
plants developed in Section 3 is robustly globally BIBO stable against not 
only unmodeled plant dynamics expressed in terms of stable factor perturba- 
tions but also uniformly bounded arbitrary internal noises as well as external 
disturbances. Although the margin of such robustness may seem pretty small 
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frora the proofs, its mere existence should be highly useful in terms of 
enhancing user's confidence on using adaptive methodologies for systems 
control in practice where a plant is seldom as imprecisely known as has been 
assumed in the proofs. 

Although the analysis of robustness margin as related to the degree of 

persistency of excitation and the adaptation gain in Section 6 is for the very 

specialized case when the regressor vector is a scalar, the result casts a 

doubt on the absolute benefit of slow adaptation. However, notice that our 

analysis did not take into account of the possible effects of varying the 

l/ 0 

adaptation gain on the uniform upper bound 0 z . True relationship between 
slow adaptation and the inherent global stability robustness margin can only 
be obtained when such effects are properly taken into account. 

It is unfortunate that, as in [4,5], the PE condition assumed in this 

paper presumes a uniform upper bound on the regressor signals before the proof 
of system stability. A recent result proving stability robustness of an AC 
methodology without a priori assuming bounded regressor is available in 
[21]. However, a complex adaptation law needs to be implemented, and 

stringent yet contradictory conditions on the adaptation gain need to be 

verified in order for the results to be valid. These conditions also favored 
slow adaptation which, as we pointed out, is not necessarily nonnegatively 
related to the inherent global stability robustness. In spite of these 

shortcomings, though, the solution methods used in [21] may still be highly 
useful for eliminating the upper-boundedness assumption made via the PE 


condition 
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